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Abstract 

The strong coupling constant qk*k-k of the scalar Kq meson decay to Kit is cal- 
culated in light cone QCD sum rule. The predicted value of the coupling constant 
9k^k-k is m a good agreement with the experimental result. 
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1 Introduction 



The number leptons that are expected to be produced yearly at planned B factories and 
the proposed r-charm factories flT|, is of the order of 10 7 , hence a detailed investigation of 
the decay properties of the r lepton becomes an important issue. The decays of r lepton 
can serve not only as a useful tool in investigation of some aspects of the standard model 
(SM) but also as a powerful experimental probe of new physics [Q. CP violation plays one 
of the most promising role in this direction. In light of this aspect, the decay of r lepton 
into hadrons has recently been investigated as probes of CP violation in the scalar sector 
of physics beyond the SM @-0. 

In [^] the Cabibbo suppressed r — ► Kirv T decay to probe the CP violation with polarized 
r's was studied. This decay mode is dominated by the contributions of the two lowest vector 
K* and scalar Kq resonances, and the mode is expected to have larger scalar contribution. 
The matrix element for r — > (Kir)~u T in the general form 

G 

M= V2 [ ^ )7m(1 ~ 75 ^ J <" + + ^)u{k)J a ] > C 1 ) 

where p and k are the r lepton and the r neutrino four momenta, respectively. The vector 
and scalar hadronic matrix elements can be parametrized as (see pj) 

Jn = sin9 c {Kit |s7 m m| 0) 



= \/2sin6 c 
J s = V2sin6 c 



F K (q 2 ) U IU - (ft - q 2 f + ^C K F s {q 2 )q, 



m 



m, — m„ 



C K F s (q 2 ) , (2) 



where 6 C is the Cabibbo angle, sin9 c = 0.23, q± and q2 are the four-momenta of 7r and K 
respectively, m s and m u are the s and u current quark masses and q = q± + q% is the four- 
momentum of the Kir system. The coupling strength Ck denoting the scalar contributions 
is determined as (see j7|) 

r _ fK* 9 Kg Kit . , 

where fx* is the leptonic decay constant of scalar Kq meson and gK*Kir is the coupling 
constant of the Kq —>■ Kn decay. In deriving Eq. (3) we assumed B(Kq — > Kir) = 100%. 
From the measured Kq Kit decay width T(Kq — > Ktt) ~ 287 MeV, the value of gK*Kn 
is obtained to be 4.87 GeV. 

In this work we employ light cone QCD sum rule to calculate Qk^k-k coupling constant 
in a model independent way and compare our results with the experimental data. 

The light cone QCD sum rule is quite different from the "classical" sum rule which is 
based on the short distance operator product expansion (OPE). This version of QCD sum 
rule is based on the OPE on the light cone, which is governed by the twist of the operators 
rather than by their dimension and the vacuum expectation values of local operators are 
replaced by the light cone hadron wave functions, and it is quite suitable in studies of the 
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hard exclusive processes in QCD. Light cone QCD sum rule has been successfully applied 
so far in the study of many different problems of hadron physics such as rare, radiative 
and semileptonic decays of B meson, S — > pj decay, nucleon magnetic moment, the strong 
couplings g n NN, Qpu-K an d 9b*b-k etc. (for an application of this method, see for example, 
the recent review || ||| and references therein). 



2 QCD sum rule for the gx^Kn coupling constant 

The aim of this section is to calculate the coupling constant gK*Kn, which characterizes the 
Kq — > Kit decay. We start by considering the two point correlation function 



U{p,q) = i / d*xe iqx (ti(p) Td(x)ij B s(x) s(0)«(0) 







(4) 



which is calculated around the light cone x 2 = 0. Here di'j^s and su are the interpolating 
currents for pseudoscalar K and the scalar Kq mesons, respectively. 

According to the basic idea of the QCD sum rule, we must calculate this correlator in 
terms of the physical particles (hadrons) and in quark-gluon language, and then equate 
both representations. 

First let us calculate the physical part of the correlator Eq. (4). Saturating this corre- 
lator by Kq and K meson states, we have 



jjphys. 



-9k*k-k 



fK S m?K S 



fKm 2 K 



(m s - m u ) (m, + m d ) f( p + q f - m 2 K ^ (q 2 - m 2 K ) 



(5) 



where (p + q) and q are the four momenta of the scalar Kq and pseudoscalar K mesons, 
respectively. In deriving the above equation we have used 







di'y^s 



K 



(K* \su\0) 

The strong coupling constant for the Kq' 

(ttK\K*) = 



fxm 2 K 
m s + m d 
fK*rn\* 

m s - m u 



(6) 



— > K°n decay is defined as follows: 

-9K*Kn ■ 



Our next task is the calculation of the theoretical part of the correlator function (4). 
The full light quark propagator with both perturbative term and contributions from vacuum 
fields can be written as 



iS(x) 



(0\T{s(x)s(0)}\0) 



i (ss) 



x- 



m 2 (ss) 



2tt 2 x 4 
9s 

167T 2 Jo 



12 192 



du< 



x 



* a afi G a \ux) - Aiu % G» u {ux) lv \ + 



(7) 
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where / = x^ . Note that here and in the following formulas the strange quark mass is 
set zero, though in numerical analysis, the mass of the strange quark is taken account. 

Substituting Eq. (7) into correlator (4) and performing Fourier transformation, for the 
theoretical part we get 

K theor - = -U Jl dtij^OOf - 4%(gi(u) + G 2 (uj) + 2g 2 (u)^ 



+ du 
Jo 



Vat; 



A? 

+2u (y\\(ai) - 2<p 



(2(p ± (ai) - (p\\(ai) + 2<p ± (aii) - <p\\(ai)) 



pq 



where 



A 
Ax 

Von 



a, 



(8) 



(m u + m d ) 

-g 2 -u - (p + qfu , 

- [q + p(ai +m« 3 )] 2 , 

dai rfa2 das S(l — a\ — a 2 — as) , 

V?(ai,a 2 ,tt3) • 



(9) 



Here ip n is the leading twist 2 distribution amplitude, ipp is the two-particle distribution 
amplitude of twist 3; g±, g 2 , <p((Xi) and <£(a:j) are the distribution amplitude of twist 4 and 



rll 

G 2 (u) = - g 2 (v)di 
Jo 



All these functions are defined as follows: 



?r(p) |di7 M 75«(0)| 0) = - f due iupx (ip{u) + x 2 gi (u)) 

J 



+ f*[ x » - ^1 / ****** 



2" > 



7r(p) ^7 5 m(0) 
t(p) d'j^ 5 g s G a(3 (ux)u(Q) 



= 



(m u + m d ) 



f 1 due iupx <p P { 
Jo 



u 



Pfl Qay. 



px 



Pa gptx 



_ %ppA 

px J 



X 



J Va tV± (a t )e ipx ^ +ua ^ , 



t(p) d~f^g s G a(3 (ux)u(0) 



Pp 5V 



XaPA ( XppA 


px J \ px J 



X 



+ i 



J Va i fi ± {a i )e ipx{ai+ua3) 

fJ}L (p aX p- PfS x a ) ( Pcw||(a;)e^ 1+ ^) , (10) 
•px J 



px 
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where the operator G a p is the dual of G a/ 3, i.e., 

G a/3 = -€ a /35 p G 5p . 

Due to the choice of the gauge x^A^(x) = 0, the path ordered gauge factor 

Vexp(ig s J dux^A^ux)) , 

has been omitted. Note that the radiative corrections to the leading twist wave functions 
are neglected, since their contribution is small (about 6-7%, see PI). 

Performing double Borel transformation with the variables (p + q) 2 and q 2 in Eqs. (5) 
and (8), we get the following sum rule for the gK*Kir coupling constant. 

9KiKirfK*fK = e f w M 

Vk* Vk 

x{ - \ <P:(u )A(s /M 2 ) + 2 g -^f (s /M 2 ) 

1 ( f u ° , F(ai,l-Mo,M -ai) 
+ M 2 \Jo 1 2(« -ai) 
f 1 da F(u , 1-u - a 3 , a 3 ) 
Jo 3 2« 3 

f u o <f\\ («!, 1 - M , m - «i) - 2(p±(ai, 1 - wo, u - Qi) 

+ / daii / r 

Jo (uq — aij 



m rl—ai \ 

da\ I — 

Juo—ai Ct 3 



(q;i, 1 - oti - a 3 , a 3 ) - 2v?±(ai, 1 - a x - a 3 , a 3 ' 



(11) 



where 



UK 
Uq 

M 2 



m, - m. 



m 



K 



m s + m d 



Mf 



(Ml + Mf ) ' 
M 2 M 2 



(Ml + Mf 



2\ ' 



with Ml and Mf are being the Borel parameters, 



dip(u) 



du 



and, 



U=UQ 



2<p ± (ai) - ip\\(ai) + 2<p±(ai) - <p\\{aii) 
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and Sq is the threshold continuum. Here the function 



n T k 
k=0 



is used to subtract the continuum and higher resonance contributions. This contribution 
is modeled by the dispersion integral, by invoking duality in the region s\, s 2 > so- in 
deriving Eq. (11), we have used the double Borel transformation formula: 



The sum rule (11) is asymmetric with respect to the Borel parameters M 2 and M| due to 
the significant mass difference of Kq and K. We choose M\ and Mf to be Mf = 2m 2 K *(5 and 
Mf = 2m 2 K f3 respectively, where (3 is a scale factor and hence in regard to this assignment 
we have M 2 = 0.44/3 GeV 2 , w = 0.107. 

3 Numerical analysis 



For numerical analysis we need the explicit forms of the wave functions. Following [|TT 
define the relevant wave functions as: 



we 



<p v (u, fi) = 6u(l - u) [l + a 2 (n)C 2 /2 (2u - 1) + a 4 (fi)Cl /2 (2u - 1)] (12) 
with the Gegenbauer polynomials 

Cl /2 {2u-l) = |[5(2n- l) 2 - 1] , 

Cl /2 (2u-l) = ^[21(2m-1) 4 -14(2u-1) 2 + 1] , (13) 
8 

and the coefficients a 2 — §, — 0.43 corresponding to the normalization point fi = 0.5 
GeV. 



<p ± (ai) = 305 2 (« 1 -« 2 )a 2 [^ + 2e(l-2« 3 )] , 

(p\\(ai) = I205 2 e(ai - a 2 )aia 2 a 3 , 

<p ± (ai) = 305 2 a 2 (l-a 3 )^ + 2£(l-2a 3 )] , 

^1 (ai) = -1206 2 a 1 a 2 a 3 [^ + e(l-3a 3 )] . (14) 



gi (u) = -5 2 u 2 u 2 + -e5 2 [uu(2 + 13uu) + 10u 3 \nu(2-3u + -u 2 ) 
2 2 5 

6 

+ 10m 3 \nu(2 - 3u + -u 2 )] , 

o 



where 5 = 0.2 GeV 2 at // = 1 Gel/ . 

The values of the main input parameteres, which appear in further numerical analysis 
are as follows: /^(l GeV) ~ 1.65, ^(1 GeV) ~ 1, / w = 132 MeF, = 156 MeV and 
m s = 155 MeV 0. 

The dependence of gx*Knfx*fK on the Borel parameter /3, for different values of the 
threshold So is given in Fig. (1). The lower bound of the Borel parameter f3 is determined by 
the requirement that, the terms of higher twists in the operator expansion must be smaller 
than the leading twist term (say 3 times). This leads to /3 > 1 for the sum rule (11). The 
upper limit of (3 is restricted from the condition that the continuum contribution must be 
less than 30% of the main one. Under this condition the upper bound is determined to be 
(3 = 1.6. 

From the analysis of Fig. (1), we finally get 



9K'K*fK*jK = 0.022 ± 0.004 



(15) 



in which we have included the uncertainties due to the continuum threshold, Borel param- 
eter, radiative corrections to the leading twist wave function, neglection of the four particle 
components of the wave functions , etc. In determination of gx*Kn the value of the leptonic 
decay constant fx* is needed. However, it should be noted that, this decay constant has 
not been measured yet, but it has been estimated in framework of different approaches, 
such as QCD sum rule which predicts fx* — 31 MeV |T3J; effective Lagrangian method 
whose estimation is ~ 45 MeV |14| and pole dominance model's result ~ 50 MeV |L5|. 
Using these values of the leptonic decay constant fx*, we get from Eq. (12) 

4.6 ± 0.8 GeV (for f K * = 31 MeV) , 



gK*Kn — \ 



3.1 ± 0.6 GeV (for f K . = 45 MeV) , 

2.8 ± 0.5 GeV (for f K . = 50 MeV) . 

Comparing these predictions with the existing experimental result gx*K-K = 4.87 GeV, it 
is observed that if we use fx* = 31 MeV, which was estimated from QCD sum rule, the 
predicted value of gx*Kn is close to the experimental result. 

In summary, we have used light cone QCD sum rule to calculate the strong coupling con- 
stant gx*K-K- The prediction we have for gx*K-K is in good agreement with the experimental 
result. 
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4 Figure captions 

Fig. 1 The dependence of gK*K-nfx*fK on the Borel parameter f3, at two fixed values of 
the continuum threshold, s = 2.4 GeV 2 and s = 2.6 GeV 2 
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